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Abstract. Observations of molecular clouds show the existence of starless, dense cores, threaded by magnetic fields. Observed 
line widths indicate these dense condensates to be embedded in a supersonically turbulent environment. Under these conditions, 
the generation of magnetic waves is inevitable. In this paper, we study the structure and support of a ID plane-parallel, self- 
gravitating slab, as a monochromatic, circularly polarized Alfven wave is injected in its central plane. Dimensional analysis 
shows that the solution must depend on three dimensionless parameters. To study the nonlinear, turbulent evolution of such a 
slab, we use ID high resolution numerical simulations. For a parameter range inspired by molecular cloud observations, we 
find the following. 1) A single source of energy injection is sufficient to force persistent supersonic turbulence over several 
hydrostatic scale heights. 2) The time averaged spatial extension of the slab is comparable to the extension of the stationary, 
analytical WKB solution. Deviations, as well as the density substructure of the slab, depend on the wave-length of the injected 
wave. 3) Energy losses are dominated by loss of Poynting-flux and increase with increasing plasma beta. 4) Good spatial 
resolution is mandatory, making similar simulations in 3D currently prohibitively expensive. 
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1. Introduction 

Magnetic fields a re observed in at lea st some molecular clouds 
JCnitcheHll999l iBourke et alJl200lh . Whether all molecular 
clouds are threaded by mag netic fields is still under debate. 
Ward -Thompson et al. ( 2000) observe ordered magnetic fields 
on small scales of about 0.05 pc in, adopting their terminol- 
ogy, prestellar cores (N ~ 10 5 cm -3 ). Also on somewhat 
larger scales , in star forming regions, ordered magnetic fi elds 
are reported JMatthews & Wilsonl2002llMatthews et alJ20 02). 
Coherent velo cities in prestellar cores are ob served on scales of 
about 0.01 pc iBarranco' & Goodmanlll998l) . On lar ger scales, 
observed line widths indicate supersonic motions. Taken to- 
gether, these observations suggest dense condensates, threaded 
by magnetic fields, to be embedded in a supersonically tur- 
bulent environment. The generation of magnetic waves under 
such conditions is inevitable. 

On larger scales, such magnetic waves are likely to be 
strongly damped (e.g. by ion-neutral friction or instabilities) 
or dominated by other processes (e. g. ISM-turbulence or in- 
coming magnetic waves as studied by Elmegreeri i 1 999l) ) . This 
finding is in agreement with molecular cloud theories and ob- 
servations. While some years ago it was thought that molecular 



Send offprint requests to: D. Folini 

* Part of this work was done during a 6 months research visit at the 
Observatoire de Strasbourg, France. 
Correspondence to: folini@astro.u-strasbg.fr 



clouds had to be supported against their self-gravity for at least 
10 8 years, new results are much more in agreement with a pic- 
ture in which molecular clouds form, stars are born, and the 
clouds are dispersed, all within some 10 6 years. Observations 
of molecular clouds i n the solar n eighborhood show that most 
clouds do form stars JHartmann et aljEoOll) . from which it is 
concluded that star formation begins essentially as soon as a 
molecular cloud forms. Using stellar evolutionary tracks leads 
to the further conclusion that star formation in a molecular 
cloud takes place rapidly, once it has started JPalla & Stahlei 
2000). For stellar populations with an average age larger than 
about 3 Myr, no mor e molecular material can be detected 
JHartmann et aljEoOll and references therein), indicating that 
star formation also ceases rapidly. N umerical sim ulations also 
support such a dynamical scenario faallesteros-Paredes et alJ 
fl999llRlme g reenl200(ilMac T,owl2002h. 

For smaller spatial scales, on the other hand, recent obser- 
vations and simulations support the idea that magnetic fields 
and waves play an important role in the structuring of the en- 
vironment of - possibly only transient - high density molec- 
ular clumps and the inhibition of accretion onto such clumps. 
Observations of th e starless dense co re L1512 and its immedi- 
ate vicinity by Falgarone et al ] j200ll) show six dense filaments 
pointing towards the core and extending up to about 1 pc. The 
matter within each filament is observed to move towards the 
core while probably describing circular motions in the direc- 
tion transverse to the filament. This motion and the orientation 
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of the filaments make it likely that they are not merely part of 
the turbulent cascade within the cloud. 

iHeitsch et alJ J200ll) performed grid studies for 3D hydro- 
dynamical and MHD simulations of the formation of persis- 
tent cores. They find that increasing spatial resolution leads to 
increased accretion in the hydrodynamical case, but to a de- 
creased one in the MHD case. The authors ascribe this differ- 
ence to better resolution of MHD waves, which then counteract 
accretion. For the case of a plane-parallel slab, at the boundary 
of which a monochromatic Alf ven wave is injected, 2D simu- 
lations with constant gravity by P runeti & Vellil dl997h. as well 
as 2D and 3D simulations without gravity by Del Z anna et alJ 
J200ll) . show the development of high-density filaments paral- 
lel to the direction of propagation of the Alfven wave. Note that 
filamentation is observed only when open, not periodic, bound- 
aries are used at t he planes perpendicular to the direction of 
wave-propagation. |Choet al. ( 2002) have shown that magnetic 
fields can have rich structures well below the viscous dissipa- 
tion scale, possibly affecting the density structure as well. 

In this paper, we study the effect of magnetic waves in 
the framework of a simplified model, a ID plane-parallel, self- 
gravitating slab with one central source of monochromatic, cir- 
cularly polarized Alfven waves. Re lated models have been in- 
vestigated by other authors before (jG arrmTie&_Ostrikerjll996t 
[ Martin et al.ll997tlFukuda & Hanawall999HFalle & Hartauistl 
l2002HKudoh & Basu 2003) 1 . The work here differs from pre- 
vious ID simulations in that we focus on the highly nonlinear, 
long-term evolution. Also, we consider only one central source 
of Alfven waves, instead of injecting energy at each grid point, 
and we use open boundaries, not periodic ones. For this setting, 
we present a dimensional analysis as well as a parameter study 
based on numerical simulations. 

Our results show that already one source of waves is suffi- 
cient to structure and support a turbulent slab in a quasi-static 
manner. The WKB solutio n for an Alfven- wave supported, 
self-gravitating ID slab by Marti n et alJJl997l) gives a good or- 
der of magnitude estimate for the average spatial extent of the 
turbulent slab, but fails to account for the rich interior structure. 
And while the Poynting-flux is constant in the WKB solution, 
there is substantial loss of Poynting-flux in the solution of the 
full equations. As governing parameters for deviations of the 
numerical solution from the analytical WKB solution we iden- 
tify the initial, central Alfven wave-length and the initial, cen- 
tral plasma beta. This is remarkable in view of the highly non- 
linear, turbulent nature of the slab, where, for example, the true 
central Alfven wave-length loses, even on average, any connec- 
tion with its initial value after a fraction of a free-fall time. To 
observe both the structuring and support of the slab by mag- 
netic waves, we find a good spatial resolution and high order of 
integration to be decisive. 

The paper is organized as follows. In Sect. |2]we describe 
our physical model and the numerical method we use. We give 
a dimensional analysis of the problem in Sect. [3] before pro- 
ceeding to the numerical results in Sect. [4] A discussion of our 
results follows in Sect.|5J conclusions are given in Sect. [6] 



2. The model 

2.1. Physical model problem 

We consider a ID (x-direction), plane-parallel, self-gravitating 
slab which we assume to be symmetric with respect to a central 
plane at x — (yz-plane, infinitely extended), where an Alfven- 
wave is injected. In this geometry, all variables are functions of 
distance x to the central plane and time t only. Velocities and 
magnetic fields perpendicular to the x-direction are allowed, 
but gradients can occur only in the x-direction. To describe the 
time evolution of this slab, we use the ideal, isothermal MHD 
equations, including a source term to account for self-gravity. 
In their conservative formulation, these equations read: 
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Here, p denotes the mass density, v = (v x 
and B = (B x , B y , B z ) the magnetic field. We will often use the 
notation B\\ for the constant background magnetic field B x , as 
well as B ± for the magnitude of the transverse magnetic field, 
and B ± o for B ± at x = and t — 0. Generally, a subscript to a 
variable refers to the initial (t - 0) and, if space dependent, cen- 
tral (x = 0) value of a variable, po is the magnetic permeability 
of free space, and /is the identity tensor. p tot = p^ + il /2po)B 2 
denotes the total pressure, where the thermal pressure is given 
by the isothermal equation of state p t h = pc 2 , with c s = y/RT 
the isothermal sound speed. T is the temperature of the gas and 
R is the gas constant. Within the framework of our ID model, 
the force exerted by self-gravity in the x-direction is given by 
g = (g(x,t),0,Q) with 



g(x,t) = -AnG p(x',t)dx'. 



Jo 



(4) 



1 The paper by Kudoh & Basu was submitted during the refereeing 
process of this paper. 



G denotes the gravitational constant. We consider only one sort 
of particle and consequently have no wave damping due to ion- 
neutral friction in our model. The average mass per particle is 
the mass mu of a hydrogen atom. 

For special cases, analytical solutions exist. For an in- 
finitely extended slab with B = 0, the stationary solution of 
the above model p roblem is given by the hydrostatic density 
distribution ( Spitze J 19681) . 

cosh (x/H) 

where po is the mass density at the central plane of the slab, 
H - c;/(27rGM hs ) is the hydrostatic scale height, and M hs = 

J-co 
„ p(x)dx is the column density of the hydrostatic slab, po and 

c s are free parameters of the solution. 

For the case where the magnetic background field By + 0, 
and when a monochromatic, circularly polarize d Alfven wave 
is injected at the central plan of the slab, Mar tin et"aD i 19971) 
derived a stationary, analytical solution of the above model 
problem in the framework of a WKB approximation. In con- 
trast to the hydrostatic solution, the WKB solution has three 
free parameters: po, c s , and B ± . 
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2.2. Model parameters and naming conventions 

The full model problem as formulated in Sect. l2.1l has five free 
parameters. They could be specified in a dimensionless form, 
as we are going to discuss in Sect. [5] Guided by observations 
of molecular clouds, we choose, however, the following set of 
dimensional parameters: the background magnetic field B\\, the 
temperature T of the slab, the initial, central mass density po, 
the amplitude of the Alfven wave, specified by either B ± o or 
v±_o = Sj_o/ V/JoPo, and the frequency o> of the wave. 

We varied parameters within limits that correspond roughly 
to observed parameters in molecular clouds: magnetic fields of 
10-100 fiG, temperatures between 5 K and 40 K, and central 
particle densities ranging from 250 to 2000cm- 3 . For the wave 
frequency u> we have assumed values in the range 10 4 yr < 
2n/o> < 2.5 ■ 10 yr. With this choice of parameters we are in a 
low-beta regime, the initial, central plasma beta Bq = 2c 2 /c 2 ^ 
lying in a range between 0.003 and 0.7. Here, cao = B\\/ sJ/ioPo 
denotes the initial, central Alfven-speed. The corresponding 
Alfven wave-length A\o = cao ■ 2jt/u>, lies in a range between 
about 0.07 pc and 0.35 pc. The magnetic field B±q corresponds 
to transverse velocities Vj_o in the range between 2.5 • 10 4 cm/s 
and 3.2 ■ 10 5 cm/s. The detailed parameters for each of the per- 
formed simulations are given in Table IA.1I of the appendix. 
They are also reflected in the name of each simulation. For 
example, R20. 10.25 .4.40 is the simulation with B\\ = 20//G, 
B ± = 10/iG, 2n/oi = 25 x 10 4 years, a central density of 
iVo = 4x 250 particles per cm 3 , and T = 40K. 

2.3. Numerical solution 

2.3.1. Numerical method 

In our simulations, we consider a slab of finite (not infinite) 
spatial extension T>. We use a finite volume method on an 
equidistant spatial grid to solve the ideal, isothermal MHD 
equations, Eqs. [QJ3] Fluxes are computed using a second or- 
der in time and third or der in space s t abilize d Lax-Friedrichs 
solver 2 as de scribed in Barmin et al. ( 1996J). As is shown in 
iBarmin et alJ Jl996l) . the accuracy of this solver is comparable 
to that of a second order Riemann solver. Self-gravity is taken 
into account using a Strang-splitting (Strang 1968). 

2.3.2. Initial conditions 

At time t — we assume the slab to have a den sity distribu- 
tion a ccording to the analytical WKB solution of M artin et alJ 
ill 9971) for a given set of parameters as specified in Table IA.fl 
The velocity in the x-direction, as well as the transverse compo- 
nents of the velocity and the magnetic field, we set to zero. The 
magnetic field in x-direction, By, and the temperature T, we set 
to the values given in Table lA.fl The mass column density of 
the initial WKB solution, within the computational domain and 
to infinity, we denote by Aio and At? respectively. 

2 T he code is part of the A-MAZE code package IWalder & Foiiml 
2000), comprising 3D adaptive mesh codes for magneto- 
hydrodynamics and radiative transfer. The codes are available at 
http://www.astro.phys.ethz.ch/stan7folini 



2.3.3. Boundary conditions 

Boundary conditions are implemented using four boundary 
cells at each of the two domain boundaries. Note that in this 
way we merely control the physical variables set in these cells. 
The fluxes entering and leaving the domain are determined by 
the interaction of the solution as set in the boundary cells and 
the numerical solution within the computational domain. 

In the boundary cells at the inner domain boundary (x = 0), 
the physical variables are set in accordance with a left-handed, 
circularly polarized Alfven-wave, whose velocity amplitude v± 
we keep fixed in time, v± = v>xo- Allowing v x and po, and thus 
B ± , to vary according to the solution in the domain, this yields: 

p(x) = p(-x), 

Vx(x) = -V x (-X), 

v y (x) = Vj_ cos(w(f - A?)), 

v z (x) - v± sin(w(f - At)), 

B x (x) = % 

By(x) = Vy(x) • ^flQp(x), 

B z (x) = v z (x) ■ y/nop(x). (6) 

For the calculation of the time retardation At = x/ca we take 
into account that ca varies with x and t. 

At the outer boundary, we distinguish two cases. If 
v x (D, t) > we use a zeroth order extrapolation for all con- 
served variables. If v x (D, t) < we use a zeroth order extrap- 
olation for v and B, but restrict the density to 10~ 4 po. Note 
that these outer boundary conditions allow for both accretion 
or loss of matter, energy, and momentum. Associated changes 
in the mass column density over the domain, Ai(t), are, how- 
ever, mostly less than 2% (see Table lA.il . 

We have chosen open boundaries as these match best with 
our intention to investigate the effect of only one source of 
Alfven waves. Using periodic boundary conditions instead 
would implicitly introduce several sources of Alfven waves, 
separated from each other by a distance 2T>. 

2.3.4. Choice of domain size, discretization, and 
integration time 

For the simulations we chose a domain of size T) — 6 pc, or 
about 20 hydrostatical scale heights, covered by 5000 cells. 
With this choice, we fulfill the following four basic require- 
ments. 1) The initial WKB solution fits well on the domain, 
90% of its column density Aio occupy less than 60% of the do- 
main. 2) The numerical solution fits well on the domain, 90% 
of its column density M{t) are contained in the inner half of the 
domain and M{t) barely changes with time. 3) The hydrostatic 
solution is covered by sufficiently many cells. 4) The numeri- 
cal solution does not depend on the discretization. Several grid 
studies show 5000 cells to be both necessary and sufficient (see 
Sect. l5.4l and Fig. [8}. We followed all our simulations for 2 • 10 7 
years, or about 20 sound crossing times of the hydrostatic scale 
height. 
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3. Dimensional analysis 

Before coming to the numerical results in Sect. |4] we present 
a dimensional analysis of the model problem formulated in 
Sect. 12.11 Dimensional parameters which completely specify 
the problem are the initial mass column density Mq, the mag- 
netic field B\\, the sound speed c s , the wave frequency a> of the 
imposed oscillations, and their velocity amplitude Vj.o at the 
central plane. Since Alfven wave propagation and self-gravity 
are essential parts of the problem, the dimensional constants no 
and G also determine the solution. The solution in infinite space 
is determined by these seven dimensional input parameters. For 
our finite computational domain this is no longer strictly true, 
but we shall neglect this complication for the moment. 

From the seven dimensional input parameters we can build 
four natural reference values for length, time, mass density, and 
magnetic field. The three extra dimensional parameters would 
then be associated with three dimensionless quantities that can 
be constructed from the seven parameters. In infinite space, the 
(suitably normalized) physical quantities in the solution would 
be functions of the normalized time and space variables and of 
these three dimensionless input parameters. 

3.1. Natural dimensional reference values 

The four natural dimensional reference values (subscript u for 
'unity' in the following) should be defined such that none of 
them would approach zero or infinity when some parameters 
of the problem take expectedly large or small values. For ex- 
ample, we should allow u> to become very large, as it might be 
in the WKB limit and as it is possibly met in actual situations. 
Similarly, c s may, in some clouds, become small enough to be 
neglected. This means that sensible natural references should 
be constructed from Mq, B\\, v±q, p-o and G alone. Obviously 
B\\ provides a reference magnetic field, B u , while the references 
of mass density, length, and time, p u , d u and f u , can be defined 
from Aio, v ± o and G alone. Straightforward dimensional anal- 
ysis shows that these four reference scales can be taken as: 



Bu = fill 



oo\2 



Pu = 



2nG(M™) 



d u = 



' ±0 



• ±0 



2nGM° 

Vj_0 

InGMl 



Pu = 



2c\ 



^2nGpu 



(8) 



The parameter a u compares the imposed velocity amplitude 
Vj_o to the reference Alfven velocity ca u , defined by c 2 Au = 
B*/(pop a ). The parameter B u is the ratio of the reference gas 
pressure to the magnetic pressure of B a - The parameter W u is a 
WKB parameter. Clearly, neither the set of the four reference 
dimensional values nor the set of the three input dimension- 
less numbers is uniquely defined. Other dimensional reference 
values could be obtained by multiplying the chosen ones by 
any function of the dimensionless input numbers a u , B u , and 
Wu- Similarly, the above set of three such numbers could be 
replaced by three other arbitrary functions of them. 

3.2. Dimensional reference values from WKB solution 

The reference quantities defined in Eq.0need not be very close 
to the actual values of the physical quantities in the solution, 
even near the central plane. It may be felt desirable, however, 
that the reference values be as close as possible to actual ones, 
even at the price of defining the reference values in a more so- 
phisticated way than p u , d u and f u . If the WKB solution is to 
be a guide, such estimates could be obtained from the solution 
of Martin et al. (1997). From this solution, it is possible to re- 
late the central density po of a solution to its column density 
Mq in the finite computational domain of thickness D and to 
the total column density to infinity, Mq . The limited extent of 
the computational domain is reflected in the fact that the ratio 
(Mq/Mq) is less than unity. A scale length d^ kb may then be 
defined as the distance from the central plane of the slab where 
90% of the column density Aio are reached. The time scale to 
can be defined as to = (2nGpoY 1 ^ 2 ■ The reference value of the 
magnetic field may be chosen as before, Bo = B\\- 

In analogy with Eq.[8] dimensionless numbers ao, Bo, and 
Wo can be introduced. Note that Bo is, in fact, identical to Bo 
as given in Sect. 12.21 An Alfven-velocity can be defined by 
c ao = B 2 /(p.oPo), and Alfven wave-length by Aaq = 2ncAo/(o. 

Aliases to the WKB parameter Wo could be used as well. In 
particular, one could prefer the ratio of the WKB length scale 
to the Alfven wavelength associated with B|| and po, Wiao = 
AaoACT ■ For our simulations, the values for a number of these 



(7) dimensionless parameters are listed in Table lA.l 



The factors 27r were inserted for convenience, such that p u and 
du are the central density and scale height of a self-gravitating 
isothermal sheet with sound speed v± that would have, in in- 
finite space, a column density Ai^. The reference time f u = 
d u /v±o is of the order of the Jeans period associated with p u . 

From the three remaining dimensional input parameters po, 
c s and to we can form, given these references, 3 natural dimen- 
sionless numbers, for example: 



(Bl/(jx oPu )) 



3.3. Reference values for substructure scale 

Unlike the WKB solution, the numerical solution of the model 
problem from Sect. l2.ll is far from smooth (see e.g. Fig. [3}- And 
while the lengths d u or c/ Q vkb are natural estimates for the global 
size of the mass distribution, the natural scale length associ- 
ated with substructures induced by the wave is, of course, quite 
different. It is either of the order of the reference Alfven wave- 
length Aau = 2jtcau/u or of the Alfven wavelength associated 
with the WKB reference quantities, Aao = 2jtcao/cl>- These two 
scales are of the same order of magnitude for a wave-supported 
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cloud. The actual characteristic scale of cloud substructure is 
the product of any one of them with a function of the dimen- 
sionless input parameters. On physical grounds, we expect, 
however, this function to be of the order of unity, since no other 
small length scale is likely to play an important role. In the low- 
P regime we consider, the sonic wavelength 2ncs /w is usually 
much shorter and sound perturbations nonlinearly evolve into 
shocks anyway, so that this sonic wavelength is not expected 
to show up in the spectrum of the fluctuations, except perhaps 
at the level of the very small thickness of the dense sheets that 
form. 

In the next section, we present the numerical solutions for 
the model problem from Sect. 12.1 1 for various parameter sets. 
Note that our study, inspired by molecular clouds, covers only 
a small part of the entire parameter space. For this part of the 
parameter space, we identify and discuss some of the depen- 
dences indicated by the dimensional analysis. 

4. Numerical results 

The injection of magnetic waves at the central plane of the self- 
gravitating slab has two major consequences: the slab becomes 
supersonically turbulent and its spatial extension is clearly 
larger than in hydrostatic equilibrium. 

Neither of these consequences is surprising. The energy 
provided by the injected wave must result in additional support 
of the slab against its self-gravity. Density inhomogeneities 
are to be expected since within the frame of linear analy- 
sis a parametric instability of the injected Alfven wave exist s 
toerbvll 19781 lGoldsteinlll978l iTurkmani & TorkelssonlEo03l) . 
Early on in our simulations we observe the associated growth 
of high density sheets, which is accompanied by the destruc- 
tion of the transver s e mag netic field (see Fig . [D. Note also 
that iMalara & Velhl dl996l) and iMalara et alTlfajQOl) demon- 
strated that even a non-monochromatic spectrum of Alfven 
waves is subject to parametric instability with linear and non- 
linear growth rates of the same order of magnitude as in the 
monochromatic case. 

So far the turbulent, nonlinear evolution of such a system 
is, however, not well investigated. The results we present in the 
following demonstrate that the same parameters which govern 
the parametric instability, A A o and fio, are also of crucial impor- 
tance for the nonlinear, turbulent solution. This despite the fact 
that the true, time averaged values of these quantities at x = 
(and close by) deviate substantially from A A o and /Jo- 

4.1. Spatial extension of turbulent slab 

We find that the time averaged spatial extension of the solution 
roughly agrees with the spatial extension of the correspond- 
ing WKB solution. As a function of time, the spatial exten- 
sion d soi (t) can oscillate but does not have to. d sol (t) we define 
in analogy with kb (see Sect. [3Jl as the distance where the 
column density reaches 90% of Ai(t). We similarly define the 
extension d hs (t) and d wkh {t) of the corresponding hydrostatical 
and WKB solutions, where corresponding means that at time t 
the three solutions have the same Ai(t). Finally, we denote by 
(.) t the time average between 5 ■ 10 6 and 2 • 10 7 years. 




0.2 0.4 0.6 0.8 1.0 

x / < d so ' > 
' t 



Fig. 1. Early on in all our simulations, high density sheets de- 
velop under the influence of the parametric instability of the 
Alfven-wave, while the transverse magnetic field is partly de- 
stroyed. The figure is a snapshot of simulation R20.20.5.4.10 
after 2.1 • 10 5 years. Shown are the density (solid line) and the 
energy density of the transverse magnetic field (dashed line), 
normalized to their initial, central values, as functions of dis- 
tance to the central plane (x-axis, in units of (rf sol )t, the time av- 
eraged slab extension). The dotted curves denote N/Nq (lower 
curve) and £'^ ag /£ , ^ iag0 (upper curve) of the initial WKB solu- 
tion. 



4.1 .1 . Dependence on system parameters 

In all our simulations, the time averaged extension (af sol ) t of the 
slab agrees to within a factor of three with the spatial extension 
of the WKB solution <f™ kb . The dominant parameters governing 
the extension of the slab are, therefore, the parameters govern- 
ing the WKB solution, i.e. po, ZJj_, and T. 

Deviations from the spatial extension of the WKB solution 
we find to depend linearly on the dimensionless WKB param- 
eter M^iao = ^Ao/<^o kb - Fig- Et shows the ratio of the spatial 
extensions <<f sol ) t /a! ( * kb as a function of A M )ld'^ h for the dif- 
ferent runs. A linear least square fit, also shown in the figure, 
gives a dependence (d sol ) t = aA M) + fid™ kh with a = 4.06 and 
/J = 0.36. Linear fits of similar quality are obtained if instead 
of c/ ( ™ kb one uses as scaling parameter <c/ hs ) t , or (c/ wkb ) t . For the 
last case, the fitting parameters are a = 3.86 and f3 = 0.37. If 
instead of A A q we consider the time average of the true Alfven 
wave-length (A A ) t at or close to the slab center, we cannot iden- 
tify any such clear dependence. (A A ) t is usually larger than A A o 
by a factor of about 2 to 8. 

From Fig. [2^ it can be seen that best agreement between 
(d sol \ and d^° is obtained around A A0 /d™ kb = 0.17. The de- 
viations at larger values of A A o/dQ kh are not too surprising. A 
basic assumption for the validity of WKB theory is that the den- 
sity changes only on scales much larger than the Alfven wave- 
length. This assumption fails to hold as A A o/d^ kb increases. For 
the deviations at small values of A A o/d^ kb we have checked that 
they are not caused by a too coarse spatial discretization, which 
would cause artificial wave-damping and thus reduced support 
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Fig. 2. a) The ratio of Aao to the extension of the WKB solution determines how well (d } t and dS agree. Shown is (d }Jd^ 
as a function of Ap&ld^ kb for the different runs. The dotted line is the corresponding linear least square fit. b) The spatial extension 
of the nonlinear solution is more or less time variable and can show distinct oscillations. For details on the oscillation period, 
see text. Shown is £/ sol (f)/c/g kb as a functions of time (x-axis, in million years) for three oscillatory runs R20.20.25.4.10 (solid), 
R100. 10.5.4. 10 (dashed), and R100.45.5.8.10 (dotted), as well as for two runs without a clear oscillation, R20.10.12.1.10 (dash- 
three-dots) and R100.45.5.4.10 (long dashes). 



against self-gravity. Increasing the spatial resolution by a fac- 
tor of four left the nonlinear solution unchanged. Instead, the 
true reason again lies in the failure of the WKB approach to 
be valid. The mass distribution in the system still consists of 
thin, dense sheets separated by more diffuse medium. As the 
sheet thickness is much less than the Alfven wavelength in the 
tenuous medium, the conditions for a WKB description are not 
met. We come back to this point in more detail in Sect. 15.11 

Besides the dominant effect of Ap^/d^ , Fig.Ek suggests 
j6o to have a second order effect. Simulations with identical ra- 
tio A^o/d^ usually have a larger extension for smaller /Jo- 
in view of what was said in Sect. [3] the average slab thick- 
ness must equal the reference scale d^ kh multiplied by a func- 
tion of ao, /3q and W.iao- In the limit of very small pressure the 
dependence on /Jo disappears. In the limit of very large Waaq, 
that is very small AM)/dQ kb , the dependence on W,\ao also dis- 
appears. There may remain some dependence on ao, which 
Fig- indicates to be weak. The ratio (d sol ) t /d^ kh then ap- 
proaches a value which is apparently close to (1/2). For higher 
values of A^o/d^ kb , we expect some dependence of (<f ol ) t /flf vkb 
on this parameter, which may be represented by a Taylor expan- 
sion to first order for not too large A^old^ . This is in rough 
agreement with what is seen in Fig. [2^, and is consistent with 
the idea that the larger the Alfven wave-length, the larger the 
spacing between dense sheets, and the thicker the system. 

4.1 .2. Time variability of slab extension 

As a function of time, the spatial extension d sol {t) of the slab 
can be very variable but does not have to be. We mention al- 
ready here that the substructure of the slab, voids and high den- 
sity sheets, always shows oscillatory motions (see Sect. l4.2l and 
Fig.|3}. In Fig. Eb, the ratio of d s0 \t)/d^ kb is shown as a func- 



tion of time for five runs, representative of the total of our sim- 
ulations. As can be seen, only some of the runs show a strong 
variability while others have a more or less constant extension. 
A measure for the strength of this time variability is the stan- 
dard deviation <x d of d &ol (t) I d^ kh . Using for the calculation of 
cr d the same time interval as for the time averages, we find the 
relative standard deviation erf = cr d /((J sol ) t /(i vkb ) to lie in an 
interval of 0.1 < erf < 0.9 (see TablelA~Tt. 

In about half of our simulations, the variability takes the 
form of a roughly periodic oscillation. The period of the os- 
cillation is, however, not always well defined. Twice the sound 
crossing time of the average slab thickness, T s = 2(d sol ) t /c s , is 
often close to the observed period r o b s . The ratio T \, s /T s aver- 
ages over the different runs to a value close to 0.95. However, 
this ratio, normalized to its average value, has significant scat- 
ter of the order of 0.8. Moreover, the fact that our slabs are more 
wave-supported than gas pressure supported makes it difficult 
to understand a scaling like r o b s ~ T s on physical grounds. A 
more natural scaling, as we will show in Sect. 15.21 would be 
that the oscillation period is proportional to twice the cross- 
ing time of the average slab thickness at velocity v_lo. For 
7vj_o = 2(£/ sol ) t /v ± o the ratio ~ T v±i) averages, over our 
different runs, to 3.4 with a scatter of 0.4. For T v±0 = 2c/ vkb /v_ L o 
the average is 3.6 with a scatter of only 0.2. Thus, to within a 
factor of order unity, but somewhat larger than unity, our simu- 
lations are also consistent with r o b s being proportional to T vl _o- 

4.2. Inhomogeneous structure of turbulent slab 

Under the influence of both the injected wave and self-gravity, 
the interior structure of the slab becomes very inhomogeneous, 
substructure develops. This substructure is the result of time- 
dependent, nonlinear effects and as such is beyond the reach of 
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Fig. 3. For simulations with the same Alfven wave-length Aao and temperature T, but otherwise widely different parameters, the 
size of the voids (black, N = 10~' cm 4 ) separating the high density sheets (white, N = 10 4 cm 4 ) is about the same. Comparison 
with Fig. |4] shows the effect of different values of Aao an d T. The general density stratification leads to larger voids at larger 
distances. Shown is the logarithmic particle density as a function of space (x-axis, in units of 6 pc) and time (y-axis, in units of 
10 6 years) for runs a) R20.20.25.4.10and d) R100.10.5.4.10. 

Fig. 4. Smaller voids (black, N = 1CT 1 cm 4 ) separating the high density sheets (white, N = 10 4 cm 4 ) are obtained for smaller 
Aao and / or higher temperatures T. a) Run R20.20.12.4.10 has the same temperature T as run R20.20.25.4.10 in Fig. [2] but 
only half the value of A A o- b) Run R20.20.25.4.40 has the same A A o as run R20.20.25.4.10 in Fig.EJbut a four times higher 
temperature T. Shown is the logarithmic particle density as a function of space (x-axis, in units of 6 pc) and time (y-axis, in units 
of 10 6 years). 



WKB theory. Again, we find Aao t0 pl av a critical role, but now 
the temperature of the slab has an effect as well. 

The density distribution is characterized by extended, low 
density voids and narrow, high density sheets. In these low den- 
sity voids, the other variables, v and B, remain nearly constant. 
Therefore, the density substructure is a good mirror for their 
substructure as well, and we restrict ourselves to density in 
the following. To characterize the density substructure in our 
gravitationally stratified slab, we represent the time series of 
ID density distributions as a 2D grey-scale plot. Low density 
regions then appear as more or less extended (dark) patches, 
while high density sheets take the form of thin (bright) lines. 

Fig- E shows such a representation of the density for two 
runs. Apparent also in this representation is the presence (run 
R100. 10.5.4. 10) or absence (run R20.20.25.4.10) of a global 
oscillation of the slab. However, it can also be seen that, inde- 
pendent of the existence of a global oscillation, the high density 
sheets within the slab undergo oscillatory motions. The motion 
of a single sheet is roughly parabolic, but may be interrupted at 
any time as two sheets collide. Different sheets can have widely 
different oscillation periods, which are mostly much smaller 
than the period of the global oscillation of the slab. 

This motion of the high density sheets, together with the 
number of sheets, determines the size of the voids (dark patches 
in Fig. |3j in the 2D representation of the data. We call this the 
scale of the substructure, a larger scale substructure thus refer- 
ring to an overall large size of the voids. Based on the only qual- 
itative measure of the grey-scale plots, we find that this scale of 
the substructure is, in essence, determined by Aao and temper- 
ature T. The size of the substructure in the two runs shown in 
Fig- El which have the same Aao and temperature T, but other- 
wise widely different parameters, is very similar. This becomes 
even more apparent if Fig. [3] is compared with Fig. 0] where 
simulations with only half the Alfven wave-length (Fig. 0Ji) 
and four times the slab temperature T (Fig.0J>) are shown. On 
the other hand, the scale of the density distribution for the two 
simulations shown in Fig.|4]is again similar. This suggest that 
augmenting the temperature by a factor of four has the same ef- 
fect as reducing the Alfven wave-length Aao by a factor of two. 
Or, put otherwise, that the scale of the density substructure in 
space and time depends on the ratio a/77 Aao or, equivalently, 
on -yySow. We currently have, however, no quantitative measure 
to further corroborate this speculation. 



The dominant role of A A o for the scale of the density distri- 
bution may have been anticipated from what was said towards 
the end of Sect.Elas well as f rom parametric instability theory 
(berbvll 19781 lGoldsteinlll97^) . Analytical results for the lin- 
ear regime indeed show the separation of density disturbances 
to increase with increasing Alfven wave-length. On the other 
hand, the system we consider is highly nonlinear and Aao soon 
loses its meaning even on average. 




Fig. 5. The energy E^ g of the transverse magnetic field is very 
inhomogeneous. It is generally smaller than in the correspond- 
ing WKB solution and not correlated with the density substruc- 
ture. Shown are the density (solid line) and the energy density 
of the transverse magnetic field (dashed line), normalized to 
their initial, central values, as functions of distance to the cen- 
tral plane (x-axis, in units (c/ sol ) t ) at the example of simulation 
R20.20.5.4.10 at a time of 2.07 • 10 7 years. The dotted curves 
denote N/No (lower curve) and ^mag/^magO ( u PP er curve) of 
the corresponding WKB solution. 



4.3. Energy of turbulent slab 

At the central plane, we constantly feed energy into the slab. 
The energy provided by this 'one point forcing' penetrates far 
into the slab, despite the gradual destruction through isother- 
mal shocks and viscous dissipation. Note that the viscous dis- 
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Fig. 6. a) Change of time averaged Poynting-flux, between the central plane and (d sol ) t , in units of the corresponding WKB 
Poynting flux, P x q as a function of -y^So for the different runs we have performed, b) Change of time averaged Poynting-flux, 
A(P x ol ) t , in units of the corresponding WKB Poynting flux, as a function of increasing, relative column density vVfoW/Alo for five 
runs, all with /3 = 0.087: R20.20.5.4.10 (dash-dotted), R20.20.10.4.10 (long dashes), R20.20.12.4.10 (dotted), R20.20.17.4.10 
(solid), R20.20.25.4.10 (dash-three-dots). 



sipation is merely given by the numerical method and the dis- 
cretization, and may not accurately mimic real viscosity, diffu- 
sion, or resistivity. The energy density of the transverse mag- 
netic field is far from smooth and is generally smaller than its 
initial WKB value. A typical situation at later times is shown 
in Fig. |5] The density substructure leads to wave reflection and 
shock formation. Occasionally, the WKB value of E^ g can be 
exceeded, as two high density sheets move against each other, 
thus compressing the transverse magnetic field between them. 

4.3.1 . Loss of Poynting-flux 

Unlike in the WKB solution, the Poynting-flux in our numer- 
ical solution is neither constant in space nor time. Let us de- 
note by P x o the Poynting-flux of the initial WKB solution and 
by (P x )t the time averaged x-component of the true Poynting- 
flux. We find that the change A(P x ol ) t of the Poynting-flux over 
a distance (d sol ) t from the central plan of the slab lies in a 
range 0.03 < A(P x ol } t /P x0 < 0.4. The ratio depends about lin- 
early on -\fjh, as can be seen in Fig. |6^. In view of the results 
in Sect. 14.21 this dependence suggests a connection between 
the loss of Poynting-flux and the amount of density structure 
we have in the slab. There we have seen that larger values of 
Vr Mao> which implies larger values of -y^So, correlate with a 
finer scale density structure. 

Largest losses of Poynting-flux occur in the innermost re- 
gion of the slab, and more than 90% of the total loss of 
Poynting-flux over the entire computational domain occurs 
within a distance {d sol ) t of the central plane of the slab. Fig. [(J) 
shows the spatial variation of the loss of Poynting-flux. There, 
(P x )t(x) - (P x ) t (x = O) is shown, in units of P x q, as a func- 
tion of increasing, relative column density Mq(x)/ Mo for five 
runs. Except for the value of a>, the runs have identical param- 
eters with a value of /Jo = 0.087. As can be seen, the decrease 
of Poynting-flux in this representation is fairly similar for all 



five simulations. Note, however, that the actual spatial exten- 
sion of the five runs is very different «t/ sol ) t = 0. 1 1 T> for 
R20.20.5.4.10but (d sol \ = 0.34 D for R20.20.25.4.10). 

The change of Poynting-flux is part of the change of the 
total energy flux. In fact, we find A(P x o/ ) t to amount to between 
60% and 90% of the change of the total energy flux between the 
central plane of the slab and (cf ol )t. The remaining amount is 
made up from the fluxes of gravitational and kinetic energy, the 
latter usually contributing considerably more. As the slab is, in 
time average, in a quasi-stationary state, this change in energy 
flux between the central plane and (flf sol ) t must be radiated. 

4.3.2. Energy equipartition 

In all our simulations, the ratio of the time averaged energy 
density of the transverse magnetic field and of the energy den- 
sity associated with the transverse velocity, (£ , ^ ag ) t /(£ , y n ) t , 
shows some spatial variation. Within about the innermost 10% 
of (d sol ) t , we find for (£ , ^ ag ) t /(£ , y n ) t typical values of about 
2, for a few simulations even somewhat higher. Further out, 
to distances around (c/ sol ) t , we observe approximate equiparti- 
tion between (E^^t and (E^ in ) t . The time averaged kinetic and 
magnetic energy densities themselves decrease with increasing 
distance from the central plane of the slab. 

If we look at the time and space averaged energy quan- 
tities, (£,t ag )t,x and (£ , y n ) tjX , we find approximate equiparti- 
tion in most of our simulations. Here, (-) t , x denotes time av- 
erage and subsequent spatial average, the later taken over 
x e [0,<£f ol > t ]. <£lg>u/<£km>t,x lies in a range between 0.8 
and 1.2 except for five runs (R100.45.5.8.10, R20.20.5.1.10, 
R20. 10.2. 1.10, R20.20.5.1.20 with a ratio between 0.6 and 0.8 
and R100.45. 1 .4. 10 with a ratio of 0.4). For the time and space 
average of the total kinetic energy density {E Km ) tyX , we find the 
ratio (£ , y n )t,x/(£kin)t,x to lie in a range between 0.7 and 0.9, 
except for run R100.20.5.4.10 where this ratio is only 0.6. 
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5. Discussion 

5.1. Spatial extension of a non-homogeneous slab 

In Sect. 0] we have seen that the WKB extension provides only 
a first approximation to the spatial extension of the fragmented 
slab. In the following, we give a qualitative explanation of why 
the extension of a fragmented slab differs from its WKB exten- 
sion, in particular also for small values of AM)/dQ kh , the seem- 
ing WKB limit. The WKB picture considers a diffuse mass dis- 
tribution. The thickness of the layer of material between col- 
umn density (the central plane) and column density M is 



H(M) 



Jr-M 




dM 



(9) 



p(M') 

where p(M) is the density of the matter sitting on top of an 
amount of material of total column density M. If a slab is frag- 
mented and contains dense sheets separated by tenuous regions, 
the contribution of the sheets to the integral in Eq. [9] is very 
small. Thus the system's thickness is essentially due to the ten- 
uous intersheet regions. This indicates that an entirely diffuse 
slab would be thicker than one which has a substantial fraction 
of its mass in the form of dense sheets, provided that the tenu- 
ous regions of the diffuse slab have drastically smaller density 
than the intersheet regions of the fragmented slab. 

In a limit where the WKB approximation is close to being 
valid, however, p(M) should be quite similar in both cases. To 
see this, let us denote by B ± (M) the magnetic amplitude of the 
wave at a point sitting on top of an amount of material of total 
column density M. Suppose the slab contains a large number of 
low-column-density thin sheets levitating in equilibrium under 
the pressure of the Alfvenic flux. Assuming equilibrium yields 
a relation between the jump dB ± of B ± at the crossing of one 
such sheet and the corresponding jump in column density dM: 

d(pA = -4nGMdM. (10) 

The same relation also holds for any infinitesimal slab of ten- 
uous material in equilibrium. Thus, irrespective of whether the 
matter has an entirely diffuse distribution or is partially frag- 
mented in a large number of thin sheets, the same relation be- 
tween B^iM) and M holds. This means that the distribution of 
wave support as a function of column density is independent of 
the details of the density distribution. If the intersheet medium 
can also be described by a WKB approach, it would result that 
the density distribution p(M) would be the same in the entirely 
tenuous slab as it is in the tenuous regions of the fragmented 
slab. A glance at Eq. [5] then shows that the fragmented slab 
would be less extended, by a factor which depends on which 
fraction of its mass is in the form of dense sheets. 

5.2. Oscillation of slab extension 



In Sect. 14.1.21 we have seen that the spatial extension of the slab 
can oscillate on a time scale of a few Myr. That a period on 
the order of T v± o = 2<t/ sol ) t /v ± o is to be expected on physical 
grounds shows the following analysis. 

We discuss the fundamental period of deviations from equi- 
librium of an isothermal slab, supported against self-gravity by 



gas pressure and a flux of WKB Alfven waves. Such a model is, 
admittedly, far from our highly structured slabs, but it should 
be sufficient to understand the main parameters which con- 
trol the global oscillation. In this model, the force exerted by 
Alfven waves takes the form of the gradient of the wave en- 
ergy density U. The wave energy flux Uca\\ we assume to be 
constant in space at any time, with cah the Alfven velocity as- 
sociated with By and the local mass density. The equilibrium is 
described by functions p e (x), U e (x) and M e (x), the latter being 
the column density between and x. An additional subscript 
denotes the values of these functions at x = 0. In particular, 
Ueo = Peov 2 0± /2. As a boundary condition for the perturbation, 
we impose, in analogy with our simulations, that the velocity 
amplitude of the wave injected at x — is fixed in time. This 
has the important consequence that the wave flux forced into 
the slab at the central plane is modulated by the varying gas 
density there. Indeed, for fixed B\\ and v_l the sum of the ki- 
netic and Poynting Alfven wave energy flux, P = pv^CA||/2, is 
proportional to ^Jp. 

In about half of our simulations, we indeed observe such a 
modulation of the Poynting-flux and density close to the cen- 
tral plane of the slab, after applying a running mean in time to 
filter out the strong short term variability stemming from indi- 
vidual high density sheets (see Fig. [7}. The modulation of the 
Poynting-flux then is in phase with the density perturbation and 
in antiphase with the slab thickness, as it should be. 

Linearizing the equations of motion about the equilibrium 
and assuming a monochromatic perturbation, M\{x)e~ ,a " of the 
column density, we find that M\{x) is a solution of the follow- 
ing linear homogeneous equation: 



p e c s 



d i 1 dM, 



dx 



dx 



U, 
+ — P, 

Pe 



3/2 



1 dMi 



dx \ n 3 / 2 dx 



, M\ (0) dU e 

2p e dx 



= 0. 



(11) 



The last term in Eg. II llrenresents the modulation in time of the 
Poynting-flux entering the slab at x = 0. It is proportional to 
the mass density perturbation at this point, M\(0) = dM/dx = 
pi(0). Had we assumed that the slab suffers perturbations un- 
der a constant wave energy flux, this term would have been 
absent from Eq. ^2 The spectrum associated with Eq. II llmav 
only be found numerically, but its gross features may be antici- 
pated from a simple dimensional analysis, substituting (ikM\) 
for (dM\/dx) and (-Uo/H e ) for (dUo/dx). H e is the thickness 
of the equilibrium slab, approximately given by the relation 



4;rGp e y/ e 2 = (t/ e o+Peo<). 



(12) 



This rough procedure does not make justice of the specific pro- 
file of p e (jt). To take care of this, we introduce numerical co- 
efficients, ct\, oc-i, q-3, q-4, all positive and of order unity. From 
Eq [n the dispersion relation is expected to be of the form: 



v 2 ikv 2 

2 2,2 , -U'0,2 . ^ lKV ±eO 

oj — ct\ c s k + CH2 k — 47r(jp e o + d\ — — . 

2 ' 2//» 



(13) 



The last, imaginary, term of Eq. ^]represents the last term of 
Eq.^2 associated with the modulation of the wave energy flux. 
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Fig. 7. The global oscillation of the slab, if existent, is in an- 
tiphase with the modulation of the running mean in time (time 
window 1.5 Myr) of the Poynting-flux and the density at the 
central plane of the slab. Shown are, as a function of time in 
10 6 years, the normalized spatial extension of the slab (dashed 
line), and the normalized running mean in time of the Poynting- 
flux (dotted line) and of the density (solid line). Normalized 
means that the quantities were divided by the maximum value 
they assume in the time interval shown. 

For the fundamental mode, k m H~ . Its frequency should then 
be given by an expression of the form: 



: 4nGp e , 



(X] 



+ (X2 



v le0 



12 



■ CUT, + ICH4 



v XeO 



(14) 



where c — c 2 + v ±e0 /2. This shows that the fundamental 
time must be, to within a factor of order unity but presumably 
smaller due to the presence of the negative term -a^, the Jeans 
frequency (AnGp e o) 1 ^ 2 . From Eq.^1 the associated time is, for 
c s « v e ox, of order of V2// e /v e o±- Given the uncertainty left 
by this analysis on the unknown factors a\ - a\, this is in rough 
agreement with the results of our simulations. The modulation 
of the Alfven flux injected in the slab, caused by the modu- 
lation of the central density perturbation, pi(0), constitutes a 
feedback which endows to 1 with an imaginary part, and thus 
induces instability, even when the system is stable under con- 
stant Alfven wave energy flux. It is interesting to note that such 
a modulation is indeed present in our simulations. The coher- 
ent oscillations observed in our calculations may represent a 
nonlinear state of development of such an instability. At high 
amplitude, the oscillations of the different fluid elements would 
cease to be isochronic and their initial coherence may be lost. 

5.3. Driving of turbulence and size of structure 

The one-point forcing we apply in our simulations is in con- 
trast to most other investigation of highly compressible tur- 
bule nce, wher e the f orcing is applied at each grid point (see 
e.g. IMacLowl i 19991) for a description of such forcing and 
Gammie & Ostriker for the particular case of a ID- 

slab). Our results show that also such one point forcing is per- 



fectly capable of supporting and structuring a slab over large 
distances. The injected waves are very efficient in distributing 
the driving energy over a large spatial range. 

For the distance over which the injected energy is spread, 
as well as for the scale of the density structure that forms, the 
driving wave-length A&o is decisive. This is not too surprising. 
In Sect. |3 we have seen that there is no other small length scale 
likely to play a role in the problem we consider. It then appears 
quite natural that A^o governs the structure size. Remembering 
that smaller scale structure means more inhomogeneity on a 
smaller spatial scale, a faster destruction of the injected energy 
for smaller Aao also seems natural. 

That local energy injection can be sufficient to drive turbu- 
lence in a whole volume has also been reported for the case 
of the interaction zone of hyper sonically colliding flows in 2D 
plane-parallel hydrodynamics JWalder & Folinilll998l l2000h . 
Although energy there is injected only at the oblique shocks 
confining the interaction zone, supersonic turbulence persists in 
the whole volume. A correlation between driving wave-length 
and structure size exists as well (Folini & Walder, in prepa- 
ration). A correlation between structures si ze and the wave- 
length of the driving was also reported by IMacLowl |[l999) 
for hydrodynamic turbulence in a 3D periodic box, which was 
monochromatically forced at each grid point. 

5.4. Numerical resolution 

We have experimented with different spatial resolutions and 
different orders of integration for the MHD equations. Not sur- 
prising, we find that simulations with shorter wave-lengths of 
the Alfven waves are more delicate in terms of spatial resolu- 
tion and order of integration. For run R20.20.5.4.10 we have 
verified that our results remain essentially unchanged (mean 
quantities differ by a few percent at most) if we further increase 
the spatial resolution. For coarser discretizations (using 1000 
cells instead of 5000) or lower orders of integration, strong 
damping of the Alfven wave due to numerical diffusion oc- 
curs and both the density substructure as well as the additional 
support against self-gravity are essentially lost, as can be seen 
in Fig. For run R20.20.5.4.10, for example, (d so1 /d™ kh \ is 
reduced from 0.46 to 0.32, for run R20.20.10.4.10 the value 
is reduced from 0.75 to 0.48. With regard to numerical simu- 
lations, these results show that the problem we investigated in 
this paper is currently out of reach for 3D simulations. 

5.5. Connection to high density cores 

We, finally, would like to return to the larger frame, molec- 
ular clouds and the formation of dense cores. As men- 
tioned earlier, both theoretical and observational results to- 
day indicate that the immediate vicinity of a forming dense 
co re is highly structured . Most striking are the observations 
bv lFalgaroneetalJ(l200ll) of the filaments around the starless 
dense core L1512. Various explanations for their existence, 
all employing magnetic fields seem plausible |Fiege & Pudritzl 
l2000t Franaueira et alJl200ol iDel Zanna et alJl200lh . It seems 
most likely that such structures, or rather the responsible phys- 
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Fig. 8. Good enough spatial resolution is crucial to observe the structuring and supporting effect of the injected Alfven-wave. 
Shown is the same simulation, R20.20.5.4.10, once, a), with our usual spatial resolution of 5000 cells and once, b), with a 
reduced resolution of only 1000 cells. Logarithmic particle density is shown as a function of space (x-axis, in units of 6 pc) and 
time (y-axis, in units of 10 6 years). A similar loss of structuring and support we observe if the order of the integration scheme is 
reduced. 



ical processes, affect the formation and accretion rate of dense 
cores. 

Detailed studies of this structuring and its consequences are 
mandatory to better link turbulence in molecular clouds with 
star formation. Such studies are, however, only at their begin- 
ning. The ID slab we studied here is by no means an accurate 
model of a molecular cloud. Nevertheless, it shows that mag- 
netic waves can play a c rucial role in c ounteracting self-gravity. 
The 3D grid studies bv lHeitsch et alJ d200lh . mentioned in the 
introduction, point in the sa me direction. And whi le the 2D 
and 3D MHD simulations bv lDel Zanna et all d200ll) of a plan 
parallel slab focus on structure formation rather than on ac- 
cretion, they show that density filaments form, more or less 
aligned with the background magnetic field. 

With regard to future simulations in this field, existing 
simulations show spatial resolution to be decisive. Another 
cruc ial issue is d imensionality. T he 'fingering' observ ed 
bv lDel Zanna etalJ feOOll) - and bv lFalgarone et alJ feOOll) - 
clearly cannot be studied in a ID model. And, as the authors 
showed as well, it can only be observed if open boundaries are 
used, not periodic ones. Finally, and probably most crucial of 
all, a future model of dense core formation will have to ac- 
count for the generation of magnetic waves self-consistently. 
This may be the strongest reason to aim, in a next step, at mul- 
tidimensional models. 

6. Conclusions 

For the case of a ID plane-parallel, self-gravitating slab we 
have shown by means of numerical simulations and dimen- 
sional analysis the following: 

1) One source of energy injection is sufficient to sustain 
turbulence throughout the slab. 

2) The time-averaged spatial extension of the turbulent slab 
is comparable to the extension of the corresponding WKB solu- 
tion. Deviations are at most a factor of three and depend about 
linearly on the ratio of the initial central Alfven wave-length 
and the extension of the WKB solution. 

3) The scale of the substructure is governed by the initial 
central Alfven wave-length and the temperature. Larger wave- 
lengths and smaller temperatures lead to larger scale structure. 

4) The energy loss, and thus the energy radiated by the slab, 
is dominated by the loss of Poynting-flux, which increases al- 
most linearly with V/6o- 

5) Within the slab, the energy density of the transverse mag- 
netic field and the kinetic energy density of the transverse ve- 
locities are in approximate equilibrium. The latter accounts for 
between 70% and 90% of the total kinetic energy density. 

(6) A too coarse mesh or a low order scheme leads to sub- 
stantial wave-damping, thus to loss of structuring and support. 
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Table A.l. The parameters used in the different simulations. The first column contains the name of each run. These names 
reflect the free parameters of each run. For example, R20. 10.25.4.40 is the simulation with B\\ = 20/j.G, B ± = lOyuG, 2tt/(jj = 
25 ■ 10 4 years, a central density of iVo = 4 ■ 250 particles per cm 3 , and T = 40K. Subsequent columns, from left to right, 
denote: A AQ = c a0 2tt/(jj: initial central Alfven wave-length in 2.176 • 10 17 cm; Alo/(ra H 2)): initial mean particle density in crrT 3 ; 
Mq/Mq'. ratio of column densities within computational domain and in infinite space of initial WKB solution; c/™ kb /c/jj s : ratio of 
extensions of WKB and hydrostatic solution at t — 0; d^ kh /D: spatial extension of initial WKB solution in units of computational 
domain; ao = B]_/B^: dimensionless parameter; fio = 2c^/c AQ : central, initial plasma beta; Wo = co/ ^2nGpo: WKB parameter; 

Waao = A A o/d™ kh : alternative WKB parameter; (d sol ),/D: time averaged spatial extension of solution, average taken between 
5 • 10 6 years and 2 • 10 7 years, in units of total domain T>\ (t/ sol ),/c/ vkb : ratio of time averaged spatial extension of solution and 
extension of corresponding, initial WKB solution; (ff: relative standard deviation of d sol (t)/d^ kb ; M20/M0: ratio of the column 
densities at times t = 20 • 10 6 years and t = 0; AP: time averaged, percental change of Poynting-flux between central plane and 
<tf ol > t in units of the initial WKB Poynting-flux, AP = 100 • (APf),/ P x0 . 
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